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ABSTRACT 

It has been recognized that magnetic reconnection process is of great impor- 
tance in high-energy astrophysics. We develop a new two-dimensional relativistic 
1 resistive magnetohydrodynamic (R 2 MHD) code, and carry out numerical simu- 

lations of magnetic reconnection. We find that outflow velocity reaches Alfven 
velocity in the inflow region, and that higher Alfven velocity provides higher 
reconnection rate. We also find Lorentz contraction plays an important role in 
O ■ enhancement of reconnection rate. 

Subject headings: MHD — magnetic fields — plasmas — relativity — pulsars: 
individual (Crab Nebula) 

c3 



1. Introduction 

Magnetic reconnection is widely recognized as a very important phenomenon in astro- 
physics. Over the last decade, it has been recognized that magnetic reconnection processes 
are very important in high-energy astrophysics. Dissipation of such super strong magnetic 
fields may play an important role both in global dynamics of the system and as a way to 
produce high-energy emission. Relativistic magnetic reconnection was proposed as a source 
of the high-energy emission (Lyubarskii 1996; Kirk et al. 2002) and as the solution to the a- 
problem (Coroniti 1990; Lyubarsky & Kirk 2001; Kirk & Skjaeraasen 2003; Lyubarsky 2003). 
Similar models were also developed for the cosmological gamma-ray bursts (Drenkhahn 2002; 
Drenkhahn & Spruit 2002; Lyutikov & Blackman 2001). Magnetic reconnection was evoked 
for explanation of the rapid variability observed in active galactic nuclei (di Matteo 1998). 
The particle acceleration in the reconnection process was proposed to operate in radio jets 
(Romanova & Lovelace 1992; Birk et al. 2001). 



- 2 - 



Due to the extreme complexity and richness of the possible effects arising in relativistic 
plasma physics, there is a strong interest for developing computer codes for relativistic mag- 
netohydrodynamics (hereafter RMHD). van Putten (1993) illustrated the implementation 
on the Riemann problem for MHD. Koide et al. (1996) then developed a RMHD code, which 
has been extensively used in relativistic two-dimensional and three-dimensional jet simula- 
tions. Komissarov (1999) and others developed and tested a Godunov-type code which is 
a truly multidimensional scheme (Balsara 2001; Koldoba et al. 2002). Recently, Del Zanna 
et al. (2003) presented a third order shock-capturing scheme for three-dimensional RMHD 
and validated it by several numerical tests. On the other hand, Koide et al. (1998, 1999) 
extended to general relativistic (GRMHD) effects, and applied it to the jet formation mech- 
anism. Gammie et al. (2003) and De Villiers & Hawley (2003) also developed GRMHD 
codes. 

Despite magnetic reconnection is recognized as an important process in high-energy 
astrophysics, there is not a lot of theoretical studies. Blackman & Field (1994) consid- 
ered kinematics of relativistic reconnection in the Sweet-Parker and Petschek configurations 
and concluded that due to the Lorentz contraction, the reconnection inflow is significantly 
enhanced and may approach the speed of light. Lyutikov & Uzdensky (2003) confirmed 
this conclusion for the Sweet-Parker case. Lyubarsky (2005) presented generalization of 
Sweet-Parker and Petschek reconnection models to the relativistic case, and argued that the 
reconnection inflow does not approach the speed of light. Particle acceleration in relativistic 
current sheets was studied both in the test particle approximation (Romanova & Lovelace 
1992; Birk et al. 2001) and in two-dimensional PIC simulations (Zenitani & Hoshino 2001; 
Jaroschek et al. 2004). Furthermore, Zenitani & Hoshino (2005) studied three-dimensional 
PIC simulations, and suggested the importance of the current-aligned magnetic field for 
studying the energetics of relativistic current sheet. Meanwhile, there are several RMHD 
simulations as we write, all these codes, however, are applied to ideal MHD and take no 
account of resistivity. 

In this paper, we develop a new two-dimensional relativistic resistive MHD (R 2 MHD) 
code, and carry out numerical simulations of two-dimensional relativistic magnetic recon- 
nection. 



2. Simulation model 

The RMHD basic equations are written as follows: 
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where c, P, v, B, E, and j are the light speed, proper gas pressure, velocity, magnetic field, 
electric field, and current density, respectively. 7 is the Lorentz factor which is defined as 
7 = [1 — (y /c) 2 ] -1 / 2 , and e is internal energy given as e = pc 2 + P/ (V — 1), where p, and T are 
the proper mass density, and the specific heat ratio, respectively. D, R, and e are defined 

as, 

D = 7 p, (6) 
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Furthermore, Ohm's law for a relativistic pair plasma under fairly general conditions has the 
MHD form (Blackman & Field 1993; Lyutikov & Uzdensky 2003), 



7 (E + " x £) = 77 J + 7 2 (j • ^ - P e c) ^ 



(9) 



where 77 is the resistivity, and p e is the electron mass density. The second term of right hand 
side of equation (9) is the convection current. 

When we solve ideal RMHD equations, it is possible to eliminate the electric field using 
E = — (v/c) x B. Therefore, past RMHD simulations solved equations (1) - (4) and evaluate 
only D, R, e, and B directly at each step from the equations. In our code, we take into 
consideration the effect of resistivity, so that we are forced to take another way. We cannot 
eliminate the electric field, therefore we also solve equation (5) to evaluate E at each time 
step. Next, from D, R, e, B, and E obtained by solving equations (1) to (5), we calculate 
7. For this purpose, we solve an equation for unknown variable 7, 



C {r( 7 2 - 1) + 1} 
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where S = (E x B)/(4irc), and P em = (P 2 + P 2 )/(87r), respectively. This equation is 
obtained by vanishing v, p and P using equations (6) to (8). We solved this equation at 
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each cell using the Newton-Raphson iteration method, so that we obtain 7. We calculate v 
and P after the iteration, and then we get p and j from equations (6) and (9), respectively. 

We assume that the evolution is two-dimensional. We take a rectangular computation 
box with two-dimensional Cartesian coordinates in the x-y plane. The medium is assumed 
to be an inviscid perfect gas. The z-component of magnetic field B z , velocity v z , and partial 
derivative d/dz are neglected. Electric field E z , and current density j z , however, are included, 
and there are no x and y components of these variables according to equations (4) and (5). 
Therefore, Ohm's law equation (9) only has z-component, and we can neglect the convection 
current term. An anomalous resistivity model is assumed, as described later. 

The region of the computation box for this study is — 52.2L < x < 52.2L, — 151.2L < 
y < 151. 2L, where L — 1 is the thickness of the initial current sheet. Non-uniform grids are 
used for both x- and y-directions. The number of grid points is 200 x 416. The minimum 
grid sizes are Ax = 0.02 and Ay = 0.05, which are concentrated near the neutral point. 

The light speed c is taken to be unity. The initial proper density outside the current 
sheet is given as p = po = 1 in non-dimensional units. We set the initial proper gas pressure 
outside the current sheet as P — Pq — 1, and the initial temperature T = P/p = Po/po = l 
is uniform everywhere. Magnitude of the magnetic field B = (8irPo/ ' (3) 1 / 2 is prescribed 
by proper gas pressure Po and the plasma (3. We study several values for (3, but it is 
P — 0.1 in the typical case. We take a relativistic Harris model as the initial current sheet 
configuration (Kirk & Skjaeraasen 2003), so we give initial conditions as B y = B tanh(2x), 
P = 1 + [13 cosh 2 (2a;)]- 1 , p=l + [(3 cosh 2 (2a;)]- 1 , v x = v y = 0, and E z = rj(dB y /dx). Since 
there is no vertical magnetic field B z , we can write E = E z z by using equations (4) and (5). 
i] is the resistivity which is defined as: 

r,( x ,y) = { ^ + ^o[2(r/r,) 3 -3(r/r T; ) 2 + l] for r < r„, (n) 

where 77^ = 5.0 x 10~ 3 is a uniform resistivity in the computation box, r] i0 = 0.3 is the 
amplitude of the anomalous resistivity, r = yj ' x 1 + y 2 is the distance from the center of the 
spot (the origin), and r v = 0.8 is the radius of the spot. 



3. Results and Discussion 

Fig. 1 shows the density distribution of the typical case ((3 = 0.1). Because of the 
enhanced resistivity around the origin, magnetic reconnection starts at this point. This 
point evolves to become an X-type neutral point. The reconnected field lines together with 
the frozen-in plasma are ejected from this X-point to the positive and negative y-directions 



- 5 - 



because of the tension force of the reconnected field lines. The velocity of the reconnection 
outflow V out — 0.9 is approximately the Alfven speed of the inflow region (Cao = 0.894). 
To complement these outflows, inflows take place from positive and negative x-directions of 
the current sheet. At the boundary between this inflow and the outflow, a shock is formed, 
emanating from the neutral point. 

Fig. 2 shows one-dimensional plots of various physical variables at t — 100, when the 
distribution becomes nearly steady state, and along y — 10, which is well upstream of the 
plasmoid ejected in the positive y-direction. At x ~ ±0.5, there are strong jumps for several 
variables. The value of current density j z becomes large, and y-component of magnetic 
field By becomes weak at these jumps. Therefore, we can say these jumps are the slow-mode 
MHD shocks. We also checked these jump conditions using the arranged model of Lyubarsky 
(2005, shown by dotted and dashed lines). From Fig. 1 and Fig. 2, we obtained tan# ~ 0.21 
where 9 is the angle between the magnetic field and the shock plane, and this value is close to 
the inflow velocity at the slow shock (e.g., Fig. 2(e)). According to the model of Lyubarsky 
(2005), inflow velocity v in ~ tan 9 in the highly relativistic regime, and our results supports 
this model. 

We next studied the dependency on the initial plasma (3. Fig. 3 shows (a) inflow 
velocity at x = 4 and y = 0, (b) maximum inflow velocity, (c) outflow velocity, and (d) 
outflow 4- velocity as function of time for j3 = 0.1, 0.2, 0.5, and 1.0 {C A o = 0.894, 0.816, 
0.667, and 0.535, respectively). Velocities are normalized by Cao in (a), (b) and (c), and 
time is normalized by Alfven transit time ta = L/Cao i n all figures. The maximum inflow 
velocity shown in (b) is almost the same as the inflow velocity at the edge of anomalous 
resistivity spot (x w ±0.8 and y — 0). Each line shows the case of a different value of (3. 
The outflow velocity reaches Cao i n all the cases. However, we obtain higher inflow velocity 
with lower (5 (higher Cao)- This means higher Cao causes higher reconnection rate v^/Cao- 
In other words, reconnection rate is higher at the relativistic regime. 

For a steady state reconnection, we can also express the reconnection rate by using the 
conservation of the mass at the steady state, V • (Dv) = V • ("fpv) = 0, so that we obtain a 
following equation: 

^ - ^Tout (12) 

v out " Pin 

where p out and p in are proper density of the inflow and the outflow region, and 7 out is 
the Lorentz factor of the outflow velocity, respectively. 5 and d are evaluated by 5/d = 
tan 9, where 9 is the angle between the y-axis and the slow-shock plane. We consider the 
Lorentz factor of the inflow velocity 7; n ~ 1. Fig. 4 shows the dependency of 5/d, p ut/Pm, 
(5/d)(p out / Pm), and 7 out to initial plasma f3 under the relativistic regime (P = 1.0), and the 
non-relativistic regime (P = 10~ 2 ). From these panels, we can see the similar behaviors of 
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ratios S/d and p ut/Pm in the both regimes. Furthermore, (S/d)(p ont / Pm) ~ 0.11 — 0.14 under 
the relativistic regime, while (S/d)(p out / p- in ) ~ 0.14 — 0.18 under the non-relativistic regime. 
Therefore the product of the two terms are roughly constant in both regimes, and we can 
say that the effect of the Lorentz factor, namely, the Lorentz contraction is an important 
factor to determine the reconnection rate under the relativistic regime. 

Let us summarize this paper. The motivation of this study is to investigate relativistic 
effects of magnetic reconnection to apply for high-energy phenomena. For this purpose, 
what we have done were; (i) to develop a new resistive relativistic MHD code, and (ii) to do 
numerical simulations for relativistic magnetic reconnection. From our study, we obtain that 
outflow velocity become close to the light speed, and due to the high inflow velocity, high 
reconnection rate is obtained. For the enhancement of the reconnection rate, we find that 
the effect of the Lorentz contraction is significant which is suggested by Blackman & Field 
(1994). However, our results also supports the suggestions of Lyubarsky (2005), specially 
in the physics of jump conditions at the slow shocks. We recognize that simulations at the 
ultra-relativistic regime are required, so we would like to report these results in future. 

The authors are greatful to Prof. T. Terasawa, Prof. M. Hoshino, Prof. S. Koide, 
Prof. S. Shibata, Dr. M. Kino and Dr. S. Zenitani for fruitful discussions. This work was 
supported by facilities of JAXA. 
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Fig. 1. — Two-dimensional density distributions of typical model (/3 = 0.1) at t = 50 and 
t = 100. Solid lines show magnetic field lines, and arrows show velocity vectors. 
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Fig. 2. — One-dimensional distributions of various physical quantities of typical model {(3 = 
0.1) parallel to the rr-axis across y — 10 at t — 100. Displayed variables are (a) proper density 
p, (b) proper gas pressure P, (c) temperature T, (d) current density j z , (e) x-component 
of velocity v x , (f) y-component of velocity v y , (g) x-component of magnetic field B x /B , 
and (h) y-component of magnetic field B y /B . Dashed lines show calculated downstream 
values, using upstream values those are shown by dotted lines based on the analytical jump 
conditions. 
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Fig. 3. — Plot of time dependence of (a) inflow velocity at (x,y) = (4,0), (b) maximum 
inflow velocity (almost same as the velocity at the edge of the anomalous resistivity spot 
(x, y) « (±0.8, 0)), (c) outflow velocity, and (d) outflow 4- velocity for several values of initial 
plasma (3. Each velocity in (a), (b) and (c) is normalized with Alfven velocity outside the 
current sheet. Time is normalized with Alfven transit time. Each line shows case for /3 — 0.1; 
solid, case for f3 = 0.2; dotted, case for (5 = 0.5; dashed and dotted, and case for (5 = 1.0; 
dashed. 
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Fig. 4. — (Left panel) Plot of plasma (3 dependence of S/d, p out / p in , (S/d)(p out / p- m ), and 7 out 
for relativistic reconnection. (Right panel) Plot of plasma (3 dependence of 5/d, Pout/Pin, 
(S/d)(p on t/ Pin), and 7 out for non-relativistic reconnection (P = 10~ 2 ). 



